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Abstract. We consider dynamo action under the combined influence of turbulence and large-scale shear in sheared 
jets. Shear can stretch turbulent magnetic field lines in such a way that even turbulent motions showing mirror 
symmetry become suitable for generation of a large-scale magnetic field. We derive the integral induction equation 
governing the behaviour of the mean field in jets. The main result is that sheared jets may generate a large scale 
magnetic field if shear is sufficiently strong. The generated mean field is mainly concentrated in a magnetic sheath 
surrounding the central region of a jet, and it exhibits sign reversals in the direction of the jet axis. Typically, 
the magnetic field in a sheath is dominated by the component along the jet that can reach equipartition with the 
kinetic energy of particles. The field in the central region of jets has a more disordered structure. 
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1. Introduction 

It is widely believed that magnetic fields play an impor- 
tant role in the formation and propagation of astrophys- 
ical jets providing an efficient mechanism of collimation 
through magnetic tension forces (see, e.g., Hughes 1991, 
Blandford 1993, Konigl & Pudritz 1999). Polarization ob- 
servations provide information on the orientation and de- 
gree of order of the magnetic field in jets. It appears 
that many jets can develop relatively highly organized 
magnetic structures (see, e.g., Cawthorne et al. 1993, 
Leppanen, Zensus & Diamond 1995, Gabuzda 1999). For 
example, radio emission from the jets of NGC 4258 in- 
dicates that the magnetic field is oriented mainly along 
the jet axis but a noticeable toroidal component can 
also be presented (Krause & Lohr 2004). The conven- 
tional estimate of the magnetic field strength in jets 
comes from the minimum energy argument and corre- 
sponds to approximate equipartition between magnetic 
and jet particle energy (see, e.g., Laing 1993). Recent ob- 
servations of jets in TeV BL Lac objects confirm this esti- 
mate (Ghiscllini, Tavecchio, & Chiaberge 2005). However, 
Kataoka & Stawarz (2005) argue that the powerful jets in 
quasars and FR II objects can be far from the minimum 
energy condition, and the field strength very likely ex- 
ceeds the equipartition value. Upper limits to the inverse 
Compton radiation of the jet in M87 imposed by HESS 
and HEGRA Cerenkov Telescopes also indicate that the 
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magnetic field cannot be weaker than the equipartition 
value (Stawarz et al. 2005). 

To explain the observational data, various simplified 
models of three-dimensional magnetic structures have 
been proposed (see, e.g., Chan & Henriksen 1980, Laing 
1981, 1999). For example, Canvin et al. (2005) computed 
emission and polarization for a few different models of 
the jets in NGC 315 and concluded that both fully and 
partially ordered fields can produce the high fractional 
polarization observed. It seems (Laing 1999) that good 
agreement with observations can be obtained if the jet 
is considered as a cylindrical core with more or less uni- 
form density and pressure surrounded by a shear layer. 
Possibly, the core and shear layer are both turbulent (see, 
e.g., Gabuzda 1999). The magnetic field has a substantial 
longitudinal component in the shear layer, but the field 
can be random or a transverse component dominates in 
the core region. Recent VLBI observations (Hirabayashi 
et al. 1998, Gabuzda, Murray, & Cronin 2004) indicate 
that the transverse field may have a substantial toroidal 
component in some objects. 

The mechanisms responsible for generation of the mag- 
netic field in jets are still unclear. Since the origin of 
jets is probably relevant to MHD-processes in magne- 
tized plasma, their magnetic fields could be generated dur- 
ing the process of jet formation (see, e.g., Blandford & 
Payne 1982, Romanova & Lovelace 1992, Koide, Shibata 
& Kudoh 1998). However, as mentioned, jets are possi- 
bly turbulent, and turbulent dissipation should lead to a 
rapid decay of any initially ordered fields if the mean-field 
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dynamo does not operate. Besides, the toroidal compo- 
nent becomes dominant as the flow expands (Begelman, 
Blandford & Rees 1984) and, under certain conditions, 
this component can be unstable thereby decollimating the 
jet (see, e.g., Spruit, Foglizzo & Stelhe 1997, Begelman 
1998). Therefore, some mechanisms for generating the 
magnetic field should certainly be operative in jets. One 
of these possibilities has been considered by Honda & 
Honda (2002), who argued that the stream can generate a 
toroidal field that participates in self-pinching the plasma 
in a fully relativistic jet consisting of electrons, positrons, 
and a small portion of ions. 

The magnetic field in jets can also be generated by 
some dynamo mechanism. It is plausible that jets are tur- 
bulent making the turbulent dynamo a suitable candidate. 
This point is argued by Stawarz et al. (2005) for the par- 
ticular case of the jet in M87 where the field strength is 
probably stronger than the equipartition value. The origin 
of turbulence can be attributed to different instabilities 
arising in jets. The classical Kelvin-Helmholtz instability 
seems to be the dominating factor of destabilization in 
the most simplified jet model with one bulk velocity and 
a narrow interface with the external medium. This insta- 
bility has been studied by many authors both in linear 
(e.g., Blandford & Pringle 1976, Ray 1981, Payne & Cohn 
1985, Zhao et al. 1992) and nonlinear regimes (Bodo et 
al. 1994, Koide, Nishikawa & Mutel 1996, Nishikava et al. 
1997, Hardee et al. 1998). In more refined models, how- 
ever, other hydrodynamic instabilities can manifest them- 
selves leading to production of turbulence in jets (Urpin 
2002a, Alloy et al. 2002). Gvaramadze et al. (1984) pro- 
posed the model where the field is generated due to the 
combined influence of helical turbulence and regular flow 
that stretches a seed magnetic field. In this model, how- 
ever, the field grows superexponcntially only during the 
initial stage but decays eventually, sothat a true dynamo 
is replaced by a temporal growth of the field. 

In this paper, we consider the turbulent dynamo ac- 
tion that can be responsible for generation of the magnetic 
field in jets. Our model is based on the shear-driven dy- 
namo action associated with turbulent shear flows (see, 
e.g., Urpin 1999, 2002b). Shear stretches turbulent mag- 
netic field lines in the direction of a mean flow, which in 
turn generates the additional component of the mean elec- 
tromotive force, which is proportional to the production of 
shear and magnetic field. Plausibly, the flow inside jets is 
sheared and turbulent, and the shear-driven dynamo can 
be in action there. There is observational evidence that the 
jet structure is rather complex and that different bulk ve- 
locities can be represented inside the jet; for more details 
see the discussion in Hanasz & Sol (1996) and Stawarz & 
Ostrowski (2002). Theoretical models of the jet formation 
also predict the existence of a nonvanishing transverse gra- 
dient in the jet velocity (see, e.g., Melia & Konigl 1989, 
Konigl & Kartje 1994, Sol, Pellctier & Asseo 1989). 

In the present paper, we show that the shear-driven dy- 
namo can naturally explain the observed magnetic struc- 
ture of jets. In Sect. 2, we derive the main equations gov- 



erning the magnetic field in sheared turbulent jets. The 
geometry and growth rate of the generated field is then 
considered in Sect. 3. A brief summary of our results is 
finally presented in Section 4. 

2. The mean-field electrodynamics of jets 

We model the jet as a cylindrical flow in the ^-direction 
with radius R. Plasma inside the jet has a velocity V = 
V(r)e z , and r, ip, z are the cylindrical coordinates with e r , 
e v , and e z the corresponding unit vectors. For the sake of 
simplicity, we assume that the jet velocity is subrelativis- 
tic. Without loss of generality, we can assume that density 
p is constant inside the jet. The Reynolds number is large 
in jets and, possibly, the plasma of sheared jets is turbu- 
lent. Variations in the flux, which are observed in many 
jets, are sometimes interpreted in terms of this turbulence 
(e.g., Marscher, Gear & Travis 1992, Massaro et al. 1999). 

We can represent the magnetic field B and the velocity 
u as a sum of the mean and fluctuating parts, B = B + b 
and u = V + v, where B and V are the mean field and 
velocity, respectively. We neglect dissipative effects in the 
induction equation and assume the field to be frozen in 
plasma. Then, averaging the induction equation, we have 
for the mean field 



ri R 

- = Vx(VxB)|Vx£, 
at 

where 
£=(vxb) 



(1) 



(2) 



is the mean electromotive force; (...) denote ensemble av- 
eraging. In what follows, however, we will need the sign of 
dissipative terms to choose the integration path properly 
when calculating Fourier integrals. Generally, a weak dissi- 
pation in jets can be provided by either electrical resistiv- 
ity (if plasma is collisional) or some plasma mechanisms, 
such as Landau damping (if plasma is collisionless). 

We consider the mean electromotive force £ using a 
quasilincar approximation. In this approximation, mean 
quantities are governed by equations including non-linear 
effects in fluctuating terms, whilst the linearized equation 
is used for the fluctuating quantities (Krause & Radlcr 
1980). A quasilincar approximation is accurate enough, 
for example, to describe an ensemble of waves with rel- 
atively high frequencies and small amplitudes. Then, the 
linearized induction equation for the fluctuating magnetic 
field reads 



— = V x (V x b) + V x (v x B). 

If V = V(r)e z , then we have 

db 8b , . 

- + V(r)--e z b r V'(r) = A, 

where 

A= (B ■ V)w — (v ■ V)B -BV-c, 



(3) 



(4) 



(5) 
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and V = dV/dr. 

Equation(4) can be solved by making use of a partial 
Fourier transformation. Since coefficients in Eq. (4) do 
not depend on ip and z we make initial transformations in 
these coordinates. The fluctuating magnetic field can be 
represented as 



b(r,t) = ^2 J dk z t>i(r,m,k z ,t)e~ 



imip—ik z z 



where m is an integer. Then, 



bi(r,m,k z ,t) 



1 



dipdze 



irrup+ik z z 



(2tt)2 

and the equation for b\ reads 

din 



b{r,t), 



at 



ik z Vbi - e z V'bi r = Ai, 



(6) 



(7) 



where Ai(r,m,k z ,t) is the corresponding Fourier ampli- 
tude of A(r, t). Substituting bi — e lkzVt b 2 , we obtain the 
equation for b 2 that does not contain the advective term, 



db 2 

dt 



^l-e z V'b 2 r = e- ik ' vt A 1 . 



(8) 



This equation can be solved by Fourier transformation in 
t. Introducing 

b(r,m,k z ,uj) = ^- j dte~ luJt b 2 (r,m, k z ,t) 

= 1 L J * e_<( " +fc " V)t &i (*"."*, ( 9 ) 
we obtain the following expression from Eq. (9) 

i - V - 

b = A(r, m,k z ,Lu) A r (r, m, k z ,u)e z . (10) 

lu lu 

Expression (10) is not a complete Fourier transform of 
b(r, t) since it depends on the radial coordinate. Note that, 
in reality, Eq. (10) does not contain singularities because 
neglected dissipative terms would result in small negative 
imaginary corrections to w, so we would have lu— iO instead 
of lu in singular terms. Here, ±«0 denotes a positive (or 
negative) contribution to lu caused by weak dissipative 
effects. The sign of this small imaginary term is important 
when calculating Fourier integrals. 

The solution for a fluctuating magnetic field reads 



dtudk 



z_^i{u+k z V)t — im<p~ik z z 



ILU 



x I A A r e z 



(11) 



Substituting this expression into the definition of £, we 
obtain 

^ J l{2TT) i U! 



e- im ^-^(v(r,t) x 



iV 

A(r u t') A r (r u t')e z 

LU 



),(12) 



where r\ — (r,ipi,z\). Taking into account that summa- 
tion over m and integration over dk z yield the correspond- 
ing (5-functions, we can simplify expression (12): 



dcudt' 
( 

2-KlLU 



,iu(t-t') 



(v(r,t)x[A(r',t') 



- — A r (r',t')e 2 

LU 



)l 



r'=r-V(t-t') 



(13) 



Integrals over dto can now be calculated using the known 
integrals if we note that lu has a small imaginary part 
caused by dissipation. We have (see Gradshtein & Ryzhik 
1965 



-f 

2TT J_, 



(-P) 



v-1 



e~ ipx dx 



{ix + P) v T(v) 



-e fjp if p < 0, 



(14) 



and if p > 0. The parameter f3 in this integral is small in 
our model because it is caused by dissipation. Then, the 
expression for £ transforms into 



£ = f dt'{v(r,t) x [A(r',t')+ 

J — CO 

(t-t')V'A r (r',t')e z }) \ r , =r _ v{t _ tl) . 



(15) 



Since A(r',t') depends on the turbulent velocity (see 
Eq. (5)) we can now calculate £ if specifying the correla- 
tion properties of turbulence. For the sake of simplicity, 
we assume turbulence to be locally isotropic and homoge- 
neous with the correlation tensor given by 

{vi(u/,kf)vj(w",k")) =Pk' i k' j S(k' + k")5(to' + to"), (16) 

where P = P(lu' } k') is the spectral function and de- 
note Cartesian components. This correlation tensor cor- 
responds to acoustic turbulence (see, e.g., Rudiger 1989), 
which seems to be plausible in a supersonic jet flow. For 
instance, sound waves can be generated by the Kelvin- 
Helmholtz instability at the jet boundary and then prop- 
agate through the jet volume generating fluctuations of 
the velocity and density (Payne & Cohn 1985). 

Since correlation tensor (16) is of particular simplic- 
ity in Cartesian components, it is convenient to represent 
the turbulent velocities v(r,t) and v(r',t') in Eq. (15) in 
terms of Fourier integrals with Cartesian wavevectors as 



v(r,t) 



dcu'dk'e 



t Auj't—ik'r 



v{k',Lu'). 



Substituting this expression into Eq. (15), we obtain, after 
ensemble averaging, 

£= f dt' f duj'dk'e^'^-^-^'^-^^iiu'.k') 

J — OO J 

x[E 1 -(t-t')V , e z x E 2 ], (17) 



where 



E 2 



E 1 = k' x [ik' 2 B + (k' ■ V')B], 
e z x k'{ik' r (k' ■ B) - ik' 2 B r - [(k' ■ V')B] r }, 



and V = {d/dr',d/rdip',d/dz'). Note that B in these 
expressions is a function of r' and t'. 
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Since turbulence is locally isotropic and homogeneous 
in a co-moving frame, the spectral function P(u)',k') 
should be an even function of the frequency u>" = 
ui' — k' z V measured in a co- moving frame, i.e. we have 
P(lu' - k' z V,k') = P(-lo' + k' z V,k'). Then, denoting 
P(lu' - k' z V,k') = G{u",k'), we can transform Eq. (17) 
into 

£= [ dLu"dk'G(Lu'\k l ) f dt'e 1 ""^ 
x [^(r', t') - (t - t')V'E 2 (r', t)]\r>=r-v{t-t>)- (18) 
Averaging of E\ and E2 over directions of k' yields 

E l = ~ V x B, E 2 = -^-e z x (VB) r , (19) 

where (VB) r = X7'B r — e ip B ip /r'. Finally, the expression 
for the mean electromotive force reads 

S = - f dt'F(t-t')[V x B 
-(t - t')V'e z x {VB) r ]\ r , =r _ v{t _ tl) , (20) 
where 

F{t) = yJ d^ dkkA G{uJ,k)e tut . (21) 

Note that the mean electromotive force given by Eq. (20) 
is nonlocal in our approach since the turbulent magnetic 
field is determined by its previous evolution under the 
influence of shear. 



3. Generation of a large-scale field in jets 

We adopt the simplest model of a jet assuming that shear 
is relatively weak in the central region but is stronger in 
a shear layer near the jet surface, r = R. If Ar is the 
thickness of a shear layer, then the shear-driven dynamo 
operates in the region R> r > R— Ar. Simulations indi- 
cate that often the thickness of a region with strong shear 
can be much smaller than the jet radius (see Alloy et al. 
1999a, b). At that point, we will not specify the radial de- 
pendence of V(r) because this dependence seems to be 
rather uncertain from both theoretical and observational 
points of view. However, we show that the proposed mech- 
anism can generate the magnetic field for any dependence 
V(r). 

In a stationary jet, the solution of the mean induction 
equation (1) can be represented as 



B{r,t) = B{r)e^- iK '*- iMv , 



(22) 



where K z is a wavevector of the magnetic field in the z- 
direction, M the azimuthal wavenumber, and 7 the growth 
rate. Solution (22) describes spiral magnetic waves. 

We are particularly interested in the generation of 
large-scale fields with not very large M. Since the thick- 
ness of a shear layer is typically smaller than R, we can 



neglect terms of the order of 1/r compared to d/dr in 
Eq. (1). Then, the r-componcnt of Eq. (1) reads 



dB r 

dt 



o r> i-oo 



/o 



(23) 



where A' is the Laplacian with the primed coordinates. 
Substituting dependence (22) for B r (r), we obtain 



/>OC 

[ 1 -iK z V{r)]B r {r)= / d^V 
Jo 



-[-y-iK,V(r)]i 



-^-KiB r + iZV'(r)K z -^ 



(24) 



where K\ = + M 2 /r 2 . Note that the radial depen- 
dence of the azimuthal component of B satisfies the same 
equation in our model. Integrating Eq. (24) over dt;, we 
have 



d 2 B r . At dB r 

~~j o - iKzV r) — 

dr z [it dr 

where 



K\ + — ) B r = 

Mr/ 



4tt _ f G(w,k) 



a + r 2 



k dojdk, 



At 



Ait 

T 



(c 2 -r 2 )G( w ,fc) ;4 



(co 2 + r 2 ) 2 



k dujdk, 



(25) 



(26) 



(27) 



and r = r(r) = 7 — iK z V(r). The coefficient /^t repre- 
sents a nonlocal magnetic viscosity in a turbulent shear 
flow, and the coefficient At describes a qualitatively new 
turbulent kinetic process that can be responsible for the 
generation of the mean field. Note that, in our nonlocal 
model, turbulent kinetic coefficients depend on the rate 
of a mean process that is the principle difference to any 
local theory like a two-scale approximation. The kinetic 
coefficients are complex in our model since T is complex. 
It is convenient to represent B r (r) as 



B r (r) - /(r)cxp (^-K z J V{r')^-dr' 



(28) 



Then, the equation for f(r) reads 

" + f -\ K ' V ' 2{r) ^ 
[it 4 [x T 



dr 2 



-K — 

2 Kz dr 



V'(r) 



At 
[i T 



f = 0. 



(29) 



The coefficients [Xt and At depend on T and, hence, 
on r, so these dependences are determined by the spectral 
function. As mentioned, the origin of turbulent motions, 
as well as their spectrum, are rather uncertain in jets. 
Turbulent motions can be caused, for example, by acoustic 
waves generated due to the Kelvin-Helmholtz instability 
(see, e.g., Payne & Cohn 1985). In this paper, we assume 
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as a certainty that turbulence is acoustic and choose the 
simplest possible dependence of G(w, k) on u), 



G(w,k) 



G{k) 

U 2 +T- 



(30) 



where t is the characteristic correlation timescalc of tur- 
bulence. Dependence (30) corresponds to a velocity corre- 
lation tensor exponentially decreasing with time, 

(vi(r,t) Vj (r,t + At)) oce- At/T (31) 
(see, e.g., Riidigcr 1989). Then, the kinetic coefficients are 



TVn 



2„,2 



At 



TV, 



1 + IV "* " (1 + Tr) 2 ' 
where vt is the characteristic turbulent velocity, 



' 2 f 



G(k)k 4 dk. 



(32) 



(33) 



Substituting expressions (32) into Eq. (29) and taking 
into account that V = —iK z V, we obtain 



dr 2 



+ « 2 (r)/ = 0, 



(34) 



where 

q (r) = K ± + - 



i K x tV" 



Ht 



4(1 + IV) 2 2 1 + IY 



To solve Eq. (34) one needs the corresponding bound- 
ary conditions. For the sake of simplicity, we assume that 
plasma inside and outside of the shear layer is highly con- 
ductive and that the mean-field does not penetrate into 
the surrounding medium. Then, function f(r) should be 
vanishing at r — R and r = R — Ar. The main qualita- 
tive conclusions are the same for other possible boundary 
conditions. 

To estimate the eigenvalues of Eq. (34), we can 
use an integral method similar to the one proposed by 
Chandrasekhar (1960). Since / is complex, we multiply 
Eq. (34) by the complex conjugate function /* and inte- 
grate over the whole region where the magnetic field is 
generated. Then, we have 



f 




2 ,-1 

dr - 


1 R—Ar 


dr 


Jr 



'R—Ar 



q 2 (r)\f\ 2 dr = 



(36) 



for the chosen boundary conditions. Both \df/dr\ 2 and 
|/| 2 are positive quantities over the integration domain, 
but q 2 (r) is complex. Splitting Eq. (36) into the real and 
imaginary parts, we obtain 



L 



R 

R—Ar 





2 „j 

dr = 


dr 


Jr 




f 




J R-Ar 



Req 2 (r)\f\ 2 dr, 



R—Ar 



Im q 2 {r) \f\ 2 dr = 0. 



(37) 
(38) 



By applying the mean value theorem, Eqs. (37) and (38) 
can be transformed into 



Re q 2 {r)dr 



R—Ar 



Ar 



Jr 



R-Ar 



Im q 2 (r)dr = 0, 



where 



Ar 


f 




(AR) 2 


J R—Ar 


dr 



dr |/(n)|- 



(39) 
(40) 

(41) 



AR <~ Ar is the characteristic radial lengthscale of the 
magnetic field, and n and r-i are some mean points within 
the shear layer. Equations (39) and (40) can be combined 
into 



q 2 (r)dr = 



R-Ar 



Ar 

TarJ 



(42) 



Substituting expression (35), we obtain the dispersion 
equation for dynamo modes 



R 
R-Ar 



3K 2 e 2 T 2 V' 2 
4(1 + IY) 2 



rv(i + rv)- 



dr 



iK z £ 2 TV" ~ 
2(1 + Tt)_ 

= -f 2 Q 2 Ar, (43) 



(35) Q 



where V" 
l 



d 2 V/dr 2 , and 



(Ai?) 5 



1 

Ar 



K\dr 



R-Ar 



(AR) 2 



M 2 
~R? 



(44) 



where Q is the characteristic wavevector of dynamo waves, 
and I = tvt is the lengthscale of turbulence. 

To estimate the eigenvalues of Eq. (43), we initially 
consider the simplest model assuming that shear is ap- 
proximately linear within the shear layer, R> r > R—Ar. 
In this case, the mean velocity can be represented as 



V(r) = V (R-r)/Ar, 



(45) 



where V is the velocity in the jet core, r < R — Ar. Then, 
Eq. (43) yields 



JrTi 



3K 2 £ 2 t 2 V' 2 



x(l + x) 



4(1 + x) 2 



dx = K z V t£ 2 Q 2 , (46) 



where T e = T(R) = 7 and T; = T(R - AR) = 7 - iK z V 
are the values of T at the outer and inner boundaries of 
the shear layer. Integrating Eq. (46), we obtain 



(1 + 7 r)( 7 - iK z V ) + % -K z V + —Q 2 - 

A T 



;K 2 Z V 2 



1 



(3£/2AR) 2 



(1+jt)(1+-/t -i T K z V ) 



0. 



(47) 



The term proportional to (3^/2Ai?) 2 is small in the 
mean-field theory and can be neglected compared to 1 
in Eq. (47). Then, the dispersion equation simplifies 



(l + 1 T)( 1 -iK z V ) + -K z V 



-K 2 V 2 + -Q 2 = 0.(48) 
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Splitting 7 into real and imaginary parts, 7 = iji + 7^, 
we obtain two equations for 7/ and 7^ from Eq. (48) 



: + , 2~„t)[ - \k x V 



0, 



(49) 



1r{1+1rt) - 7/r(7/ - K Z V ) - T -K 2 Z V 2 + — Q 2 = 0.(50) 

6 T 



Equation (49) yields 

11 = \k z v q . 

Then, we have from Eq. (50) 

(T7«) 2 +T 7 *-(^X-W 
The roots of this equation are 



0. 




— K 2 V 2 



(51) 



(52) 



(53) 



One of the roots is positive, so the corresponding dynamo 
mode is growing if 



-K 2 V 2 > £ 2 Q 2 . 



(54) 



T 

12 

Estimating Q 2 w l/(Ai?) 2 + K 2 and assuming that the 
azimuthal wavelength is much larger than AR, we can 
represent Eq. (54) as 

1/2 

(55) 



V 



K Z AR — > 2V3 ( 1 

VT 



\2v\ 



This condition cannot be fulfilled if Vo < 2\/3vt, and 
the dynamo does not operate in such jets. However, it is 
plausible that the mean velocity is much larger than the 
turbulent velocity in jets, Vo 3> 2^/3vt, and condition (55) 
can be satisfied for a wide range of dynamo modes with 



(56) 



where X z = 2tt/K z is the wavelength in the ^-direction. 
Therefore, the maximum longitudinal wavelength of the 
dynamo-generated magnetic structure in jets is of the or- 
der of AR(V /v T ). 

The generation time, t* = I/7.R, is given by 

-1 



i* = 2t 




(V 2 124) 



U 2 

(KkJ 



.(57) 



For dynamo waves with a relatively short wavelength, 
\ z < 2tt£(Vq/vt) (or Vq ^> 1/tK z ), the growth time is 



2V3 



(58) 



K z V a 

In the limit of a large wavelength, X z > 2tt£(Vo/vt) (but 
still satisfying condition (55)), Eq. (57) yields 

U ~ T7T^ TTo- (59) 



Kit 2 V 2 ' 



The dynamo modes with the maximum possible wave- 
length (~ AR(Vq/vt)) grow on the slowest timescale 
U ~ t(AR/£) 2 . 

It is seen from this consideration that generation of the 
mean-field is determined by a velocity difference between 
the boundaries of a shear layer rather than by the details 
of the velocity profile. Therefore, our results can be gener- 
alized for any velocity profile in a relatively simple way. We 
can split Eq. (43) into the real and imaginary parts and 
then obtain equations analogous to Eqs. (51) and (52) by 
applying the mean value theorem. The only difference to 
Eqs. (51) and (52) is that the equations for a more com- 
plicated velocity profile will contain the value of a flow 
velocity on some mean point V(r*) (R > r* > R — Ar) 
instead of Vo- Correspondingly, the growth rate and the 
generation condition will be given by Eqs. (53) and (54), 
respectively, with the replacement Vo — > V(r*). Therefore, 
a particular shape of the velocity profile appears to be 
unimportant in our model, and large-scale magnetic fields 
can be generated in any jets with a strong shear. 



4. Discussion 

We have considered the turbulent dynamo action in jets. 
The main result is that even the simplest turbulent mo- 
tions showing the mirror symmetry become suitable for 
the generation of a large-scale magnetic field in the pres- 
ence of shear. An amplification of the mean field takes 
place due to non-local terms that appear in the mean 
electromotive force and are caused by shear stresses. The 
considered mechanism of generation is qualitatively differ- 
ent from the conventional turbulent alpha-dynamo that, 
apart from the lack of the mirror symmetry of turbulence, 
also requires large-scale stratification. Due to its simplic- 
ity, the proposed mechanism is well adopted to the phys- 
ical conditions in jets because the presence of both shear 
and turbulence seems to be plausible in a jet flow. 

Unfortunately, neither available observational data nor 
theoretical modelling provide reliable information con- 
cerning the velocity profile in jets. However, the generation 
of a large-scale magnetic field can take place for any ve- 
locity profile, which is an attractive feature of our model. 
The only necessary condition of the considered dynamo is 
the presence of a sufficiently strong shear satisfying condi- 
tion (54). This condition can be fulfilled in many jets or, 
at least, in a fraction of their volume. This dynamo mech- 
anism generates the field in jets on a very short timescale 
that can be comparable to the turnover time of turbu- 
lence r and that is typically much shorter than the lifetime 
of jets. Therefore, the generation can most likely reach a 
saturation level when the dynamo works in the nonlinear 
regime. 

We considered generation of the radial component of 
B since this component is the most important one in the 
shear-driven dynamo. Generally, two other components 
can be stronger but their evolution is determined entirely 
by the behaviour of B r . We can estimate B z from the z- 
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component of the mean induction equation, equating the 
dissipative and stretching terms. Then, 

B z ~ V '( AK)2 B r ~ K Z AR ( — } B r , (60) 
Mt \vtJ 

and the longitudinal magnetic field is always stronger 
than the radial one. The azimuthal magnetic field can 
be estimated from the divergence condition and can vary 
within a wide range depending on the longitudinal wave- 
length and M. For not very large B v is typically 
stronger than the radial field. Our linear analysis does not 
allow proper estimation of the saturation magnetic field 
but, most likely, the strongest field component can reach 
cquipartition with the kinetic energy of jet particles. 

In our model, a dynamo generates a large-scale field in 
the layer of the thickness Ar near the jet surface. The 
shear-driven dynamo is much less efficient in the core 
region, r < R — Ar. However, if the jet is turbulent, 
the small-scale turbulent dynamo can amplify small-scale 
magnetic fields in the core region even if the shear is neg- 
ligible. The magnetic Reynolds number is large in jets, 
and turbulent motions caused, for example, by instabili- 
ties stretch and distort the field lines increasing thereby 
the energy of generated small-scale magnetic fields rapidly 
(see, e.g., Schekochihin et al. 2001). This random field 
reaches equipartition with the energy of turbulent motions 
on a very short timescale <~ r. Since turbulent motions 
are less energetic than the mean flow, we can expect that 
random fields in the central region are typically weaker 
than the large-scale magnetic field in a magnetic sheath 
surrounding the jet. The characteristic lengthscalc of tur- 
bulent magnetic fields in the central region can be shorter 
than i. 
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